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ABSTRACT 

The dispersion in lithium abundance at fixed effective temperature in young cool 
stars like the Pleiades has proved a difficult challenge for stellar evolution theory. We 
propose that Li abundances relative to a mean temperature trend, rather than the 
absolute abundances, should be used to analyze the spread in abundance. We present 
evidence that the dispersion in Li equivalent widths at fixed color in cool single Pleiades 
stars is at least partially caused by stellar atmosphere effects (most likely departures 
from ionization predictions of model photospheres) rather than being completely 
explained by genuine abundance differences. We find that effective temperature 
estimates from different colors yield systematically different values for active stars. 
There is also a significant correlation between stellar activity and Li excess, but not a 
one-to-one mapping between unprojected stellar rotation (from photometric periods) 
and Li excess. Thus, it is unlikely that rotation is the main cause for the dispersion 
in the Li abundances. Finally, there is a signficant correlation between detrended Li 
excess and potassium excess but not calcium- perhaps supporting incomplete radiative 
transfer calculations (and overionization effects in particular) as an important source 
of the Li scatter. Other mechanisms, such as very small metallicity variations and 



^Current address 



- 2 - 



magnetic fields, which influence PMS Li burning may also play a role. Finally, we 
find no statistical evidence for a decrease in dispersion in the coolest Pleiades stars, 
contrary to some previous work. 



Subject headings: open clusters and associations: individual (Pleiades) — stars: 
abundances, activity, atmospheres, interiors, late-type, rotation 



Introduction 



Predictions are made by standard stellar models (e.g., Bahcall Sz Ulrich 1988 ) about the 
surface abundances of elements in stars. However, there are indications that such models are 
incomplete. A case in point is the surface abundance of the element lithium (Li) in low mass stars, 
which is observed to decrease with time. 



The solar meteoritic value for the Li abundance is 3.31ib0.04 (Anders & Grevesse 1989). 
A study of Li abundances in young, pre-main-sequence (PMS) T Tauri Stars (TTS) suggests a 
value of log N(Li) = 3.2±0.3 ( |Magazzu, Rebolo, &: Pavlenko 19921) , consistent with the meteoritic 
value. While this is one indicator of the initial Li abundance, TTS abundance determinations are 
beset by complications due to their youth, such as uncertain Tgfj estimates and the presence of a 
circumstellar accretion disk. The study of stars in open clusters of different ages, like a Persei (50 
Myr) and the Pleiades (70-100 Myr), shows that there is nearly a uniform Li abundance of 3.2 for 
high mass stars (~7000 K). 

It is known that surface Li depletion takes place during the PMS evolution of low mass stars 
due to Li burning via {p, a) reactions at low temperatures of T;^2.6xl0^ K. Surface depletion 
can occur in standard models through convective mixing if the base of the convection zone is 



hot enough to burn Li (Bodenheimer 1965, Pinsonneault 1997). Because PMS stars have deep 
convection zones, they burn Li during the PMS. As the depth of the convection zone is a function 
of mass (increasing with lower mass) Li is depleted on the main sequence only in lower mass 
stars (<0.9 Mq). However, the Pleiades evinces a large dispersion in surface Li abundance at a 
given color for Tcfj^5500 K (e.g., Soderblom et al. 1993b| ). Standard stellar models are unable 
to reproduce this dispersion. Furthermore, open cluster observations indicate some depletion is 
observed on the main sequence as well, which is in conflict with standard standard models. 

Because these models cannot fully explain the observed depletion patterns, additional mixing 
mechanisms seem necessary. Rotation provides one driving mechanism for such non-convective 
mixing, through meridional circulation ( [Tassoul 1978 , Zahn 1992 ) and instabilities caused by 



differential rotation (Zahn 1983). Hence, rotation in stars has received much scrutiny as a possible 
agent of Li depletion and of the observed scatter in open cluster Li abundances at a given mass. 
Models which include rotational mixing ( Pinsonneault, Kawaler, Sz Demarque 1990| ) are able 
to predict the dispersion seen in older systems, but not at young ages like that of the Pleiades 
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(Chaboyer, Demarque, & Pinsonneault 1995). The study of Li abundances is a rich and vast field, 
and there have been several efforts to study the correlation of surface Li abundances with rotation 
using stars in open clusters. Here, we concentrate on the connection between surface Li abundance 
and rotation using data in the young Pleiades cluster. 

Because the Pleiades Li scatter is such a difficult obstacle in our understanding of early stellar 
evolution, a historical summary seems in order. Butler et al. (1987) studied a sample of 11 K-stars 
in the Pleiades and determined that four rapid rotators had higher Li abundances than four slow 
rotators. They believed this consistent with the evolutionary picture that on arrival on the main 
sequence, stars had high rotation rates and high Li abundances (i.e., they arrived on the main 
sequence before there was time for rotational braking or Li depletion). As the star spun down, the 
Li abundance decreased as well. Hence, they concluded that the faster rotators were younger than 
the slower, hence less depleted. 

A study of the distribution of rotational velocities of low-mass stars in the Pleiades by Stauffer 
& Hartmann (1987) revealed that there was a wide range of rotation velocities in the Pleiades K 
and M dwarfs. They showed that the distribution of rotation velocities in the Pleiades could be 
reproduced quite well invoking angular momentum loss, without having to resort to a large age 
spread which is also in conflict with the narrow main-sequence seen among the low-mass Pleiades 
stars. 

Soderblom et al. (1993b) carried out an extensive study of Li abundances in the Pleiades. 
They considered several possible explanations for the dispersion in the observed abundances, 
including observational errors and the effect of starspots. They concluded that the spread in Li 
abundances seen was real and not an artifact of other physical conditions. They found that the 
Li abundance was correlated well with both rotation and chromospheric activity, and speculated 
that rapid rotation was somehow able to preserve Li in stars. While they found some low vsini 
systems with high Li abundances, it was possible that these stars are faster rotators simply seen 
at low inclination. 

Balachandran et al. (1988) studied a sample of stars in the younger a Persei cluster (50 
Myr) and concurred with the picture of Li-poor stars as slow rotators. However, a comparison 
of the Pleiades to a Per by Soderblom et al. (1993b) showed that while most stars had similar 
abundances, a significant number of stars in a Per had abundances that were less than that in 
the Pleiades by 1 dex or more. This was difficult to understand until Balachandran et al. (1996) 
published a corrected list of Li abundances that culled all non-members from the sample, bringing 
consistency to the Pleiades and a Per abundances. 

Garcia Lopez et al. ( 1991a, b) added seven stars to the Butler et al. sample in the range 
4500;^Tofj;^5500 K and asserted a clear connection between Li abundance and vsini. They also 
found that the correlation breaks down for temperatures cooler than 4500 K. In a subsequent 
paper in 1994, they enlarged their sample and further studied the correlation, concluding that 
their earlier assertion was correct - there were no rapid rotators with low Li abundances and there 
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was a clear relationship between log N(Li) and vsini. They did note three stars (H II 320, 380 
and 1124) having low vsini values and Li abundances comparable to those of the rapid rotators 
as counterexamples, but speculated these objects were rapid rotators seen at low inclination 
angles. It is to be noted that, when determining mean abundances for their rapid and slow rotator 
populations, they included the 3 stars with low vsini and high Li abundances in their rapid 
rotator sample. While this does not change the qualitative result they obtained, it does affect the 
magnitude of the difference in abundance between the two populations. It also demonstrates that 
there is a range of abundances for slow rotators. 

Jones et al. (1996a) derived Li abundances and rotation velocities for 15 late-K Pleiades 
dwarfs, and also found that the correlation between Li abundance and rapid rotation breaks 
down for cooler stars (Tefj^4400 K). Jones et al. (1997) determined rotational velocities and Li 
abundances in the 250 Myr old cluster M34, intermediate in age between the young Pleiades 
(70-100 Myr) and the Hyades (500 Myr). They concluded that the Li depletion and rotation 
velocities were in between the Pleiades and Hyades values, and that the pattern seen in these 
clusters suggested an evolutionary sequence for angular momentum loss and Li depletion. 

One could speculate that some of the Li dispersion in the Pleiades and a Per may be due 
to NLTE effects and unknown effects of stellar activity on the Li I line formation ( [Houdebine &: 



Doyle 1995 ; Russell 1996| ; but see Soderblom et al. 1993b). The structural effects of rotation might 



also be responsible for the Li depletion pattern. Martin & Claret (1996) included this ingredient 
in their models for masses of 0.7 and 0.8 M© and were able to produce enhanced Li abundances 
for stars with high initial angular momenta as a result of less effective PMS Li destruction in 
rapid rotators (i.e., their models imply initially rapid rotators will have high Li abundances 
relative to the other stars at the same T^ff at young ages). Angular momentum loss as well as 
rotationally-induced mixing could affect these models. 

The questions we consider here are if there truly is a correlation between Li abundance and 
rotation rate in the Pleiades, what the nature of the correlation is, and if not, what might explain 
the abundance scatter. We start with a careful sample selection for our analysis, and examine 
various possible causes that might contribute to the dispersion (including errors in abundance 
determination). We then proceed to an analysis of the Li-rotation correlation and explore other 
possible correlations that might be masquerading as a Li-rotation correlation. 



2. Pleiades Lithium Abundances 

2.1. Sample Selection and Definition 

Lithium abundances were derived from the datasets in the studies of Soderblom et al. (1993a; 
S93a) and Jones et al. (1996). The two studies were merged to form our starting sample with 
the latter data preferred in cases of overlap. Secondary stellar companions can affect photometric 



-5- 



colors from which T^e values are derived, activity levels, measured line strengths, and activity 
levels deduced spectroscopically. Thus, in order to look at the intrinsic Pleiades Li abundance 
dispersion unrelated (directly or indirectly) to the presence of a stellar companion, binary systems 
were excised from our sample. Cluster and interloping field binaries identified by Mermilliod et 
al. (1992) and Bouvier et al. (1997) were removed from the starting sample. Two spectroscopic 
binaries not in these lists, but identified as such by S93a, were also removed. 

The color-magnitude diagram of this refined sample was then inspected to photometrically 
identify binaries using the dereddened BVI photometry described by Pinsonneault et al. (1998). 
We found H II 739 to be an obviously overluminous (or overly red) outliar in the V ys. B — V,V 
vs. V — I, and / vs. V — I diagrams, and eliminated it from the refined sample. Finally, all stars 
with upper limits on the A6707 Li I line's equivalent width were eliminated. These upper limits, 
as censored data, complicate the ensuing statistical analysis. These stars are also the very hottest 
and very coolest in the sample. Their photometrically-inferred Teff values and model atmospheres 
may be slightly more uncertain than the other objects in the sample. Their elimination simplifies 
the analysis and reduces possible additional sources of uncertainty. 

This final sample of 76 Pleiades stars is listed in the first column of Table 1. The 

extinction-corrected V magnitude and reddening-corrected {B — V) and iy — I) colors are given 

in the second, third, and fourth columns. The color-magnitude diagram of these stars evinces a 

tight main sequence, and is shown in Figure 1 (open circles) with a 100 Myr isochrone described 

in Pinsonneault et al. (1998) and assuming a distance modulus of 5.63 (Pinsonneault et al. 1998). 

The only possibly discrepant outliars remaining are: a) H II 686, which appears overluminous in 

the \ YS. B — V plane, but not the V — I plane b) H II 676, which appears underluminous (or too 

blue) in the V — I plane and perhaps B — V also, and c) H II 2034, which appears underluminous 

(or too blue) in the B — V plane, but not m V ~ I. There is no convincing evidence that these 

slight discrepancies are related to binarity. Rather, they may be due to relatively large photometric 

errors in one passband or to physical effects {e.g., increased red flux from spots) unrelated to 

binarity. Stars rejected as binaries axe plotted as filled triangles; their general propensity to reside 

above the main sequence is evident. Tab. 1 

Fig. 1 

2.2. Stellar Tgs and Activity Measures 

S93a and Jones et al. (1996) provide photometric T^s estimates for all of the stars in Table 
1. We re-examine these for comparison and because of concern that chromospheric activity or 
starspots might affect the colors of young stars. We calculated T^s using our {B — V)o values and 
the relation from Soderblom et al. (1993b, equation 3): T^s = 1808(S -y)^- 6103(5 -y)o + 8899. 
Temperatures were also derived from our (V — I)o colors using the relation from Randich et al. 
(1997): Teff = 9900 - 8598(y - 7)o + 4246(F - 1)1 - 755{V - I)^. Both of these relations are 
based on the data from Bessell (1979), and should provide self-consistent temperatures given 
self-consistent photospheric colors. Columns 5, 6 and 9 give the Tgs values of S93a, and those 
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derived here from (B — V) and (V — I). 

We adopt the Ha- and Ca II infrared triplet-based chromospheric emission measurements 
from S93a as stellar activity indicators. These are the ratio of the flux (relative to an inactive 
star of similar color) in the Ha and Ca II lines relative to the total stellar bolometric flux. 
Given canonical views of a relation between stellar mass and chromospheric emission on the 
main sequence, it is also of interest to measure the residual Ha and Ca II flux ratios. That 
is, we wish to detrend the general relation between stellar mass and activity such that activity 
differences unrelated to large-scale mass differences can be quantified. This was done by fitting 
the Ha and Ca II flux ratios as a function of {V — I) color temperature with a linear relation^, 
and subtracting this fitted flux ratio (computed at a given V — I) from the measured flux 
ratio of each star. The relation for the fitted Ha flux ratio (used below) was found to be 
logi2Ha,fit = (-0.00044742xTefr(y - I)) - 2.12515. The relation for the Ca II flux ratio was found 
to be logi?caii,fit = (-0.00021017xreff(F - I)) - 3.50280. 

We find strong evidence that our T^g values (hence, assuming self-consistency of the color- Tgfj 
relations, the photometric colors) are affected by activity level. Figure 2 shows the difference 
between the {B — V)- and (V — /)-based T^s values versus the Ha flux ratios (top panel) and the 
mass-independent residual Ha flux ratios (bottom panel). A relation is seen in both panels, such 
that the lowest T^g differences are seen predominantly for the lowest flux ratios while the largest 
Teff differences are seen predominantly for stars having the largest flux ratios. The ordinary linear 
correlation coefficients are significant above the 99.9% confidence levels for both panels. 

The binary stars (filled triangles) behave similarly to the single stars in both panels; on 
average, though, the binaries exhibit larger Teg residuals than the single stars. This systematic 
offset likely reflects the additional influence of fainter (hence cooler and redder) companions on the 
photometric colors. If this interpretation is correct, it could suggest that the handful of inactive 
single stars with significant T^s residuals in the upper left portion of both panels are unrecognized 
binaries. 

Significant differences between the {B — V)- and (V — I)-based Tgg estimates, the slight 
propensity for (B — V) to yield larger temperatures, and the association of these properties with 
stellar activity seems to be a common property of young stars noted and discussed by others {e.g., 
Garcia Lopez et al. 1994; Randich et al. 1997; King 1998; Soderblom et al. 1999). Explanations 
for these observed properties in young stars are at least twofold: a) increased -B-band flux due 
to boundary layer emission associated with a circumstellar disk (presumably not applicable for 
our near-ZAMS Pleiads), and b) increased I-band flux due to the presence of cool spots. It is 
straightforward to associate increased prevalence and surface coverage of spots with increasing 
activity, and Ha emission (used here to quantify activity) has been associated with accretion of 



^Quantitative comparison of the resulting values indicated that fits with higher order functions did not yield 
statistically improved descriptions of the flux ratio-color relations. 
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circumstellar material in young stars. Given that the Pleiades age (~100 Myr) is an order of 
magnitude larger than inferred disk lifetimes for solar- type stars (Skrutskie et al. 1990), spots are 
the more likely cause of the temperature difference-activity relation in our Pleiades sample. 

In their recent study of the effects of activity on Pleiades Li abundances, Stuik et al. (1997) find 
that activity- specifically the presence of spots and plages- may significantly alter photospheric 
colors. Indeed, they suggest that the resulting changes in color may be a more dominant 
contributor to the Pleiades Li spread than line strength differences. Additionally, they find 
that such color variations are both surprising and complex. Their empirical solar-based activity 
models indicated that both spots and plages lead to increased {B — V) colors; in contrast, their 
"best-effort" theoretical stellar models indicate a decrease in [B — V). Sorting out which (if either) 
set of models are appropriate for specific Pleiads (in addition to other empirical details such as 
specific spot/plage coverage and ratio) might be further illuminating. Because the direction of 
changes in {V — I) also flips in their models, Stuik et al. (1997) note that spot/plage-related 
changes in color may not be ideally identified in two-color plots (e.g., Soderblom et al. 1993). 



2.3. Lithium Abundance Determinations and Detrending 

Li abundances for all our Pleiads were determined from the measured A6707 line strengths 
and our preferred Tgfj value. A blending complex lies some 0.4 A blueward of the Li doublet. In 
our stars, the typical contribution of these blending features is '^10 mA, which is significantly 
smaller than typical Li line strength of ~100 mA. The blending contribution was subtracted^ 
following the approach of S93a, who parameterized the contaminating line strength as a function 
of {B — V) color. Here, we recast this parameterization as a function of Tgg so that differences in 
our {B — V)- and {V — /)-based temperatures were consistently accounted for in the analysis. 

Given the Tefj values and the corrected Li line strengths, abundances were calculated using 
Table 2 from S93a. This was done by fitting a surface map of the equivalent width-temperature- 
abundance grid of S93a using high order polynomials. The internal interpolation accuracy 
is generally a few thousandths of a dexQ. Columns 8 and 11 of Table 1 give the derived Li 



^Deblending corrections were not applied to any stars taken from Jones et al. (1996) following their claim that 
instrumental resolution was sufficient to separate the Li line and blending complex. While it is not clear to us that 
this is true given that some of their objects have appreciable rotation (see their Figure 1), it does not affect the present 
results inasmuch as the Jones et al. rapid rotators have Li line strengths significantly larger than those expected of 
the blending complex. 

^The two hottest single stars have T^g values significantly outside the curve of growth grids provided by S93a. 
Extrapolation to these temperatures with high order polynomials leads to errantly low abundances by a few tenths 
of a dex. To the extent that we are mainly interested in the cooler Pleiads and that we are only interested in the 
differential star-to-star Li abundances (i.e., large scale abundance morphology with T^s is removed later), these known 
errors are unimportant for the present analysis. In the case of these two stars (H II 133 and 470), we simply caution 
those who would use our absolute abundances, and also note that the few much smaller extrapolations to lower Tcff 
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abundance^ for our {B — y)-based T^g, and for our {V — /)-based T^g. 

At the Pleiades age, PMS Li burning has significantly depleted the initial photospheric Li 
content of many of our stars. Moreover, this PMS depletion is a sensitive function of mass (or 
Teff) with less massive stars having depleted more Li due to deeper convection zones and longer 
PMS evolutionary timescales. In examining star-to-star Li abundance differences connected with 
parameters such as rotation or activity, we must remove this general large-scale trend in the 
abundance vs. T^s plane. 

The procedure is illustrated in Figures 3 and 4, which shows the Pleiades Li abundances 
versus our {B — V)- and {V — I)-based Tq^. The familiar and large (3 orders of magnitude) 
abundance depletion over 2000 K of T^g is seen in Figure 3. The mean trend is shown by the 
dashed line, which is a fourth order Legendre polynomial fitted to the single star data after 
rejecting ±3a outliars. Fourth order fits were also conducted for the data based on the S93a 
temperatures and our {B — y)-based values. These fits provide a mean fiducial Li abundance at 
any given T^g to which observed abundances (calculated assuming the same source of T^g) can be 
compared to infer and measure a relative Li "enhancement" or "depletion" for each star as shown 
in Figure 4. For {V — /)-based temperatures, the approximation to the fit shown in Figure 3 is 
given by: log iV(Li)fit = -10.8602 + (2.9785x 10"^ xT) + (1.1736x lO^'^xT^) - (3.8181x 10"" xT^). 
For (B — y)-based temperatures, the approximation to the fit shown in Figure 3 is given by: 
logAf(Li)fit = -10.9959 + (2.8360x 10"^ xT) + (1.7354x lO'^xT^) - (4.2744x 10~" xT^) 

2.4. Errors 

Uncertainties in the Li abundances were estimated from those in T^g and equivalent width. 
Here, we are only interested in the internal errors which affect the star-to-star Li abundances. A 
measure of the internal uncertainties in the T^g estimates is provided by the estimates from the 
(B — V) and (V — I) colors. For the single stars, the difference in the two color-temperatures 
exhibits a per star standard deviation of 108 K. Assuming equal contributions from both {B — V) 
and {V — I), this suggests an internal error of ±76 K in the T^g of any one Pleiad derived from 
any one color. This uncertainty was translated to a Li abundance error by re-deriving abundances 
with Teg departures of this size. The adoption of identical errors in (B-V) and (V-I)-based Teff 
values is a simplifying assumption (though one likely true within a few tens of K). Inasmuch as 
our conclusions are the same using either the (B-V) or (V-I) colors, it is not a critical one for this 
work. 

The other significant source of uncertainty is in the Li line measurements. For S93a line 
strengths (the majority of our sample), uncertainties come from their own quality measures: a 



outside the S93a grid are not believed to be affected by any substantial amount, 
^by number, relative to hydrogen, on the usual scale with log A*'(H)= 12. 
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(±12 mA), b (±18 mA), c (±25 mA), and d (±40 mA). Jones et al. (1996) state that their 
uncertainties range from 5-20 mA and depend largely on projected rotational velocity. Assuming 
this range and the stars' vsini values, we have adopted the reasonable values shown in column 12 
of Table 1. Given S93a's note that the equivalent widths of Butler et al. (1987) may have to be 
regarded with caution, and the typical expected uncertainties from Poisson noise expected from 
their S/N and resolution, we have assigned an uncertainty of ±30 mA in these line strengths. 
Based on the S/N, resolution, and vsini values of the observations from Boesgaard et al. (1988), 
we have adopted a conservative uncertainty of ±5 mA for their data. In a similar fashion, we 
assigned uncertainties of ±25 mA to the equivalent widths from Pilachowski et al. (1987). The 
line strength uncertainties were translated to Li uncertainties by re-deriving abundances with the 
adopted equivalent width departures. 

Final Li abundance uncertainties, shown in Figures 3 and 4, are calculated by summing the 
two errors in quadrature, and listed in columns 13 and 14 of Table 1. We emphasize that for 
the purpose of looking at the star-to-star Li abundance differences in cool Pleiads, the effects of 
Teff errors are minimized. This is because the movement of a star in the T^s-Li plane due to T^s 
errors is very nearly along the cool star depletion trend for Tcfj^5800 K. To take into account 
this correlation in looking at the differential Li abundances (i.e., the actual values versus an 
expected value from a fitted trend to the data), the abundance errors due to departures in T^g 
were combined with the slope of the fitted Li versus T^s trend at the Tgs of each star. The total 
uncertainties in the differential Li abundances are given in the final two columns of Table 1. 



2.4.1. Li Abundance Scatter 

Large scatter in the star-to-star Li abundances is apparent in Figures 3 and 4. Comparison 
of the observed scatter with that expected from the estimated uncertainties indicates that the 
spread is statistically significant. The presence of real global scatter was considered by comparing 
the variance of the differential Li abundances with that based on the uncertainties given in Table 
1. The sizable reduced chi-squared statistic (Xu — 12.78, v = 72) indicates probabilities of the 
observed variance (s(Li)^~0.13 dex^) occurring by chance are infinitesimal. 

Additional analysis was carried out by binning in T^g. For both the {B — V) and {V — I) 

based results, we broke the data up into 5 Tefj ranges following natural breaks in the estimated 
Tcff values which yielded comparable sample sizes (10-15 stars) in each bin. The results for both 
the (B — V) and (V — I) data are similar. 

We find that stars in the hottest bin (bin 'E': 6172-6984 K and 6107-6928 K ior B - V and 
V — I) exhibit a variance that is larger than the expected value at only the 94% confidence level. 
The stars in the adjacent cooler bin (bin 'D': 5567-6048 K and 5521-6021 K for B - F and F - /) 
exhibit a variance in the Li abundances significantly larger than expected from the uncertainties 
at the ;^99.93% confidence level. The differential Li abundances in the remaining three cooler bins 
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(bin 'C: 4899-5477 K and 4996-5452 K; bin 'B': 4507-4815 K and 4542-4746 K; bin 'A': 3955-4332 
K and 3867-4343 K) all show observed variances significant at considerably higher confidence 
levels. 

An important claim by Jones et al. (1996) in their study of Pleiades Li abundances is a 
progressive decline in the dispersion of the Li abundances as one proceeds from the late G dwarfs, 
to the early-to-mid K dwarfs, and finally to the later K dwarfs. We find, however, that quantitative 
analysis fails to provide firm support for such a conclusion. F-tests of the observed variances 
indicate that differences of the differential Li abundance dispersions in our cooler three bins are 
statistically indistinguishable. The differences between the bin B and bin A stars' variances are 
significant at only the 71.5% and 78.0% confidence levels for the B — V and V — I datasets. 
Differences between the bin C and bin B stars' variances are different at only the 72.7% and 
75.0% confidence levels. It should be noted that these comparisons ignore the observed Li upper 
limits prevalent for the coolest (bin A) Pleiads. The stars with upper limits lie at the lower edge 
of the observed Li abundances (figure 4 of Jones et al). Ignoring this censored data may lead 
to an underestimate of the true dispersion for the coolest Pleiads- maJiing our conclusion of no 
significant difference in the magnitude of star-to-star abundance scatter for the late G to late 
K Pleiads a conservative one. Larger samples and improved upper limits (or detections) would 
clarify this important issue. 



3. Nature of the Li-rotation correlation 

3.1. Projected rotational velocity 

Extant studies of Li-rotation correlations have employed v sin i measurements, which yield 
only a lower limit to the rotational velocity due to the unknown angle of inclination, i. In Figure 
5, our V — /-based absolute and differential Li abundances are plotted against the projected 
velocity measurement t;sini. The top two panels (a and b) show the data for all Teff. The bottom 
two panels (c and d) show data with 4500<Tcff <5500 K, which is the temperature range in which 
a clear connection between Li abundance and vsini was asserted by Garcia Lopez et al. (1994). 
It is seen that while there is a range of abundances at the lower values of the rotation velocity, 
the rapid rotators (t^sinz;^ 30 km s~^) do show a tendency to have higher Li abundances in the 
intermediate T^e range (panel d). 



3.2. Rotational Period 

The Pleiades now has many members with photometrically-determined rotation periods 
(Krishnamurthi et al. 1998), which are free of the ambiguity associated with inclination angle. 
Hence, it is now possible to consider the true nature of the correlation between rotation and Li 
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abundance. For example, we find a slowly rotating star (H II 263, P = 4.8 d) that has a high 

Li abundance. Furthermore, two of the three stars (H II 320 and 1124) in the Garcia Lopez et 

al. (1994) study with low vsvui and high Li abundances also have measured rotation periods of 

4.58d and 6.05d respectively. Thus, there are several cases where high Li abundance in stars with 

low vsm.i is not due to inclination angle effects- Li overabundances are not solely restricted to 

rapid rotators. This is apparent in Figure 6, where the surface Li abundance is plotted against 

rotation period, Prot; instead of vsvui. In particular, we draw attention to the large range in 

abundances seen at longer rotation periods (>4.0 days; panels b and d). Thus there exist at least 

a few Pleiads which are true slow rotators, but have high Li abundances. Fig. 6 

We next examined the proposal by Garcia Lopez et al. (1994) that there is a very clear 
relationship between rotation and log N(Li) for stars with M ~0.7-0.9 Mq. Figure 7 shows the 
V — /-based Li abundances versus mass for Pleiads with photometrically-measured rotation rates. 
The symbol size is proportional to the rotation period. When rotational periods are considered 
rather than v sin i measurements, a true range of Li abundances with rotation is seen in the mass 
range 0.7-0.9 M©- there are genuine slow rotators with abundances similar to the fast rotators. 
Thus, there appears to be a range of rotation at all abundances. Hence, Prot is essential to study 
the true correlation. Fig. 7 



3.3. Structural effects of rotation 

Rapid rotation affects the structure of a star and hence the derived mass at a given Tefj 
( Endal &: Sofia 1979| ). The structural effects of rotation would alter a star's color such that rapidly 
rotating objects would be redder, hence perceived as cooler, and thus be assigned a lower mass. 
We therefore examined the abundances as a function of mass rather than effective temperature. 

To investigate this issue, it was necessary to construct stellar models for different disk 
lifetimes ( |Krishnamurthi et al. 1997| ). We ran models for uJcrit = ^xloq and Wxujq to represent 
fast rotators and the slow rotators. The rotation-corrected masses were derived by interpolation 
in the models across effective temperature and rotation velocity for different disk lifetimes. We 
found that the percent change in mass is small (5%) even for the most rapidly rotating star in the 
Pleiades (H II 1883, usini=140 km s^^). The change is between 1% and 2% for stars with vs'mi 
in the 50-100 km s^^, and < 1% for vsini < 50 km s~^. These small alterations fail to eliminate 
the Li dispersion, which sets in at M<0.9 Mq in the Li-mass plane. Thus, the structural effects of 
rotation on the derived mass-temperature relation are not large enough to account for the Pleiades 
Li abundance dispersion. 

These results differ with those of Martin & Claret (1996), who also explored the structural 
effects of rotation and found enhanced Li abundances for stars with high initial angular momentum. 
This is not seen in our models, which predict small rotational structure effects on the Pleiades Li 
abundances, similar to the models of Pinsonneault et al. (1990). Mendes et al. (1999) have noted 
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the conflict between the results of Martin & Claret (1996) and Pinsonneault et al. (1990), and 
considered the hydrostatic effects of rotation on stellar structure and Li depletion using their own 
stellar models. Their results are in agreement with Pinsonneault et al. (1990), and they too find 
that hydrostatic effects are too small to explain the observed Li abundance spread in the Pleiades. 



4. Li and Stellar Activity 

4.1. Li and chromospheric emission 

Since the large Pleiades Li spread is in such a young cluster, one may wonder if its 
long-sought explanation is related to stellar activity. Additionally, since rotation and activity are 
well-correlated, a Li-activity relation may be masquerading as a Li-rotation relation instead. Here, 
we discuss if magnetic activity indicators such as chromospheric emission (CE) are correlated 
with the Li abundance. Several studies have pointed out that activity is correlated with the Li 
abundance (e.g., ^oderblom et al. 1993b , Jones, Fischer fc Stauffer 1996b| ). There have also been 



some studies speculating that CE affects the apparent abundance of Li (e.g., [Boudebine fc Doyle 



1995 ) and hence may be at least partly responsible for the dispersion. 



Figure 8 contains our results, and shows the V — /-based differential Li abundances versus 
the Ca II infrared triplet fluxes (top panel) and residual fluxes (bottom panel). A relation is seen 
in both panels, such that the lowest log N(Li) differences are seen predominantly for the lowest 
flux ratios while the largest log N(Li) differences are seen predominantly for stars having the 
largest flux ratios. The ordinary correlation coefficients are significant at the 99.7% and >99.9% 
confidence levels for the chromospheric Ca fluxes and residual fluxes, suggesting a significant 
relation between chromospheric activity differences and Li abundance differences (though not 
necessarily causal). 



4.2. Spreads in Other Elements 

Important clues to the cause of the Pleiades Li abundance scatter can be found from 
examination of other elements not destroyed in stellar interiors like ^Li. Variations in such 
abundances may signal effects other than differential Li processing, and perhaps point to an 
illusory difference caused by inadequate treatment of line formation. 



4-2.1. Potassium 

One of the most useful features for this purpose is the A7699 K I line. The usefulness of this 
feature is two-fold. First, there is the similarity in electronic configuration with the Li I atom and 
the fact that this particular K transition and the A6707 Li I are both neutral resonance features. 
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Second, the interplay of abundance and ionization effects leads to the happy circumstance that 
the line strengths of these two features are comparable in Pleiades dwarfs. Thus, line formation 
for both features should be similar in many respects. 

The K I line strengths were taken from S93a and Jones et al. (1996). These were 
then plotted versus Tcff as derived from both B ~ V and V — I. The relation was 
well fit by a 4th order Legendre polynomial approximated by the relation: EW(K) = 
9098.151 - (4.126458XT) + (6.381173x IQ-^xT^) - (3.308942x IQ-^xT^) for the B - F colors and 
by the relation: EW(K) = 8926.171-(4.170605xr) + (6.702972xl0-^xr2)-(3.642286xl0-8xr3) 
for the V — I colors. These fits showed considerable scatter - the line strength dispersion was ~55 
mA, which is considerably larger (and statistically significant) than even the maximum equivalent 
width errors estimated by S93a. So scatter is present in the potassium data as well. 

Differential K I equivalent widths ( [observed— fitted] /fitted ) are plotted against the 
differential Li abundances in Figure 9. A correlation between the values is present, though with 
considerable scatter. The one-sided correlation coefficients are significant at the 99.0 and 98.3% 
confidence levels for the {B — V) and {V — /)-based results. Such a correlation (of whatever 
magnitude), however, may arise not from some physical mechanism; instead, it may simply reflect 
correlated measurement errors. Fig. 9 

Like the differential Li abundances, the differential K I line strengths are correlated with 
activity measure. Figure 10 shows the differential K I equivalent width versus the Ca II fluxes (top 
panel) and residual fluxes (bottom panel). The correlations arc analogous to those seen for the 
differential Li abundances in Figure 8, and are significant at the ~98.5% (top panel) and >99.9% 
(bottom panel) confidence levels. Fig. 10 



4.2.2. Calcium 

To examine the possibility of correlated measurement errors, we considered the 
line strengths of the A6717 Ca I feature taken from S93a. These were fitted against 
Teff in the same manner as the K I equivalent widths. The relations are given by: 
EW(Ca) = 4203.706- (1.859218xr) + (2.888087x10-4 xT^) - (1.538325xl0-8xr3) for theB-V 
colors and by EW(Ca) = 2978.740 - (1.249348 x T) + (1.887209x lO'^xT^) - (9.973237x lO^^xT^) 
for the V — I colors. The scatter associated with these fits is ~20 mA, which is consistent with the 
S93a uncertainties. Interestingly, unlike Li and K, there is no evidence for scatter in the calcium 
data above the measurement uncertainties. Differential Ca I equivalent widths are plotted against 
the differential Li abundances in Figure 11. The relation is flat. Unlike K, there is no significant 
correlation- the ordinary correlation coefficients are significant at only the ~80% confidence level. Fig. 11 



This indicates to us that the correlated scatter in Li and K line strengths is not due to 
measurement errors. Rather, we suggest that some physical mechanism affecting the details of 
line formation not included in standard LTE model photosphere analyses is the cause. Such a 
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mechanism, if having star-to-star differences, may be the dominant source of the Li abundance 
scatter in Pleiades dwarfs. Since activity evinces such differences, it may naturahy provide such a 
mechanism. 

In an important theoretical study, Stuik et al. (1997) have urged similar caution in regarding 
Pleiades Li scatter as solely due to genuine abundance differences. These authors consider the 
photospheric effects of activity on Pleiades Li I and K I line strengths by modeling surface spots 
and plages. They can neither exclude nor confirm these particular manifestations of magnetic 
activity as the cause of the problematic and important K I variations in cool Pleiads. Their 
extensive efforts, though, do open the door for future improvement. 

First, it seems important to establish whether their empirical solar-based spot/plage models 
or their "best-effort" theoretical stellar models are more nearly correct, and if one or the other 
model set is indeed applicable to all Pleiads since the two model sets produce color and line 
strength changes opposite in sign. Second, Stuik et al. (1997) note that their radiative equilibrium 
and mixed activity calculations depart from observations with increasing [B — V). As they 
acknowledge, such disparities may signal other effects not yet considered: UV "line haze", which 
may depend on the presence and structure of an overlying chromosphere, impacting the details of 
line formation; unknown properties and effects of Pleiads' granulation patterns; and the influence 
of so-called solar-like "abnormal granulation" within plages. Third, other sources of non-thermal 
heating of the photosphere by chromospheric 'activity' may need to be considered. Finally, simply 
relating colors or effective temperatures (from color-Tgg conversions) of Pleaids having different 
activity levels may be more problematic than realized. 

Houdebine & Doyle (1995; HD95) demonstrate that formation of the A6707 Li I line is 
sensitive to activity in M dwarfs. The extent of these effects depends on the relative coverage of 
plages and spots. HD95 note the particular importance of the role of ionization in reducing the 
resonance line's optical depth. In late G and K dwarfs like those showing scatter in the Pleiades, 
star-to-star variations in departures of both photoionisation and collisional ionization from that 
predicted by model photospheres might introduce significant star-to-star variations in the derived 
Li abundance^. Interestingly, King et al. (1999) find element-to-element abundance differences in 
two cool (Teff~4500 K) Pleiades dwarfs and a similarly cool NGC 2264 PMS member which are 
ionization potential dependent. We suggest that current evidence may implicate non-photospheric 
ionization differences as a likely source of star-to-star Li variations in the Pleiades. 



5. Other Mechanisms and Concerns 



Overionization from photospheric convective inhomogeneities has been discussed in the context of Population II 
star Li abundances by Kurucz (1995). 



- 15 - 



5.1. Metal abundance variations 

Variations in abundances of other elements can affect stellar Li depletion via the effects 
of stellar structure on PMS Li burning. For example, Figure 3 of Chaboyer, Demarque, & 
Pinsonneault (1995), indicates that very small metal abundance differences of, say, 0.03 dex lead 
to substantial (;^0.3 - 0.4 dex) differences in PMS Li burning for Tcfr^4500 K. 



Extant studies ( |Boesgaard fc Friel 1990 ; Cayrel, Cayrel, &: Campbell 1988 ) of Pleiades F- and 



G-star iron abundances (which cannot simply be equated with "metallicity" when it comes to 



PMS Li depletion; Swenson et al. 1994) suggest no intrinsic scatter larger than 0.06-0.10 dex. The 
photometric scatter of the single stars in the color-magnitude diagram might allow a metallicity 
(or, perhaps more properly, those elements which are dominant electron donors in the stellar 
photospheres) spread of 0.05 dex. While small, these constraints would still permit substantial Li 
abundance spreads for cool Pleiads. Additionally, abundances of elements which may have a large 
impact on PMS Li depletion but little effect on atmospheric opacity (e.g., oxygen) have yet to be 
determined in cool Pleiades dwarfs. 

However, the Li spread in the Pleiades extends to Tefj values substantially hotter than ~4500 
K. At hotter Tgff values, model PMS Li-burning is less sensitive to metallicity. For example, in 
the range 5000-5200 K, the observed Li abundance spread would require "metallicity" differences 
approaching a factor of two. Such spreads would be surprising indeed, and not expected based 
on the limited results of extant Fe analyses of hotter cluster dwarfs. Abundance differences (of a 
large number of elements) of this size would not be difficult to exclude with good quality spectra 
as part of future studies. 



5.2. Magnetic Fields 



Ventura et al. (1998) have recently investigated the effects of magnetic fields in stellar models 
and PMS Li depletion. They find that even small fields are able to inhibit convection, and thus 
PMS Li depletion. They suggest that a dynamo generated magnetic field linked to rotational 
velocity (thus, presumably yielding an association between activity and rotation given conventional 
wisdom) would result in ZAMS star-to-star Li variations that mirror differences in star-to-star 
rotational (and presumably activity) differences. As these authors admit, the fits of their magnetic 
models to the Li-Tegp morphology and significant scatter of the Pleiades observations are not 
"perfect" or "definitive"; however, the qualitative agreement and ability to produce star-to-star 
scatter and general relations between Li abundance and rotation and activity are encouraging. 
Continued observations (especially detailed spectroscopic abundance determinations of various 
elements in numerous Pleiads) and theoretical work will be needed to establish the degree to which 
the Pleiades Li spread is illusory or real and, if the latter, its cause(s). 



6. Summary and Conclusions 



The very large dispersion in Li abundances at fixed T^g in cool (r;^5400 K) Pleiads is a 
fundamental challenge for stellar evolution because standard stellar models of uniform age and 
abundance are unable to reproduce it. A variety of mechanisms (rotation, activity, magnetic fields, 
and incomplete knowledge of line formation) have been proposed to account for this scatter. Here, 
we construct a sample of likely single Pleiads and consider this problem with: a) differential Li 
abundances relative to a mean Tes trend b) rotational periods instead of projected rotational 
velocities c) chromospheric emission indicators, and d) line strengths of other elements. 

We calculated Tgs values from both {B — V) and {V — I) on a self-consistent scale based 
on the calibrations of Bessell (1979). We find differences in the two T^s values and these are 
significantly correlated with both general activity level and with differences in activity, suggesting 
that surface inhomogeneities may noticeably affect stellar colors. Our results are consistent with 
a growing body of evidence of significant differences between (B — V)- and (V — /)-based Teff 
values, a propensity for {B — V) to yield larger T^s values, and a relation of these characteristics 
with activity in young stars from 5 Myr old PMS stars in NGC 2264 (Soderblom et al. 1999) to 
PMS stars in the 30 Myr old IC 2602 (Randich et al. 1997) to ZAMS stars in the ~100 Myr old 
Pleiades. 

However, the similarity between the sensitivity of the derived Li abundance to TeS and the 
clusters' physical Li-Tefj morphology means that even substantial Tcfj errors are not a significant 
source of star-to-star Li scatter. Nor are observational errors. Comparison of the scatter in 
the differential Li abundances with errors from T^s and line strength uncertainties indicates an 
infinitesimal probability that the observed scatter occurs by chance. We find significant scatter in 
the Li abundances below ~6000 K; it is significantly larger, though, below ~5500 K. Statistical 
analysis fails to support previous claims of smaller scatter in the late K dwarfs relative to the late 
G and early-mid K dwarfs. 

There is a spread of Li abundance at low vsini, whereas the rapid projected rotators tend to 

have larger differential Li abundances in the range 4500<Teff <5500. However, use of photometric 
rotation periods (free from uncertainties in the inclination angle i) indicates there is not a 
one-to-one mapping between differential Li abundance and rotation. The stars H II 263, 320, and 
1124 are examples of stars with Li excesses but slow rotation (P= 4.6 — 6.1 d). In contrast to 
previous claims based on vsini, the rotation periods indicate a true range of Li abundance with 
rotation in the mass bin 0.7-0.9 M0. 

Using the theoretical framework of Krishnamurthi et al. (1997), we constructed stellar 
models to investigate the hydrostatic effects of rotation on stellar structure and PMS Li burning. 
As shown in Figure 8, these models fail to account for the Pleiades Li dispersion, which is in 
agreement with the independent conclusions of Mendes et al. (1999). 

We find that the star-to-star differences in Pleiades Li abundances are correlated with activity 
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differences, as measured from Ca II infrared triplet flux ratios, at a statistically significant level. 

Moreover, the Li differences are significantly correlated with differences in the strengths of the 
A7699 K I resonance feature. This seems to not be due to correlated measurement errors since the 
Li differences show no correlation with the A6717 Ca I line strength residuals. This is a significant 
result given similarities in the Li and K feature's atomic properties and line strengths. We suggest 
that incomplete treatment of line formation, related to activity differences, plays a significant role 
in the Li dispersion- i.e., that part of the dispersion is illusory. As emphasized by Houdebine & 
Doyle (1995), the formation of the Li I feature is sensitive to ionization conditions. If chromosphcric 
activity variations can produce significant variations in photo- and collisional-ionization in the Li I 
line formation region not accounted for by LTE analyses using model photospheres, this may lead 
to errors in the inferred abundance. Relatedly, we note the results of King et al. (1999) who found 
ionization potential-dependent effects in the elemental abundances of two cool (Tes) Pleiads and a 
similarly cool NGC 2264 PMS star. 

If such conjecture is correct, we expect that somewhat older (less active) cluster stars will 
exhibit less Li dispersion. This seems to be the case for the ~800 Myr old Hyades cluster 
(Thorburn et al. 1993) and perhaps also for the mid-G to mid-K stars in M34 (Jones et al. 1997). 
These clusters still exhibit scatter, and this may be real and due to differences in depletion from 
structural effects of rotation (Mendes et al. 1999), magnetic fields (Ventura et al. 1998), small 
metallicity variations (§5.1), main sequence depletion due to angular momentum transport from 
spin-down (Pinsonneault et al. 1990; Charbonnel et al. 1992) or a planetary companion (Cochran 
ct al. 1997), and photospheric accretion of circumstellar/planctary material (Alexander 1967; 
Gonzalez 1998). The amount of scatter expected in even older clusters is less clear. If, e.g., 
rotationally-induccd mixing acts over longer timescales then the scatter may well increase again; 
indeed, the substantial Li scatter in M67 solar-type stars observed by Jones et al. (1999) could 
indicate that this is the case. 

These authors called attention to the possible pattern of very large Li scatter in young 
clusters, considerably reduced scatter in intermediate-age clusters, and increased scatter in older 
clusters. In the scenario we envision, variations in activity-regulated ionization of the Li I atom 
may be responsible for the majority of (mostly illusory) star-to-star Li differences in near-ZAMS 
and younger stars; of course, this does not exclude a (lesser) role from other variable mechanisms 
influencing PMS Li burning. If the decline in the activity level of intermediate-age stars reduces 
the importance of variable ionization, then the (smaller) Li scatter in these stars could arise 
from variations in PMS Li burning due to, e.g., the hydrostatic effects of rotation on stellar 
structure, inhibition of convection by magnetic fields, and small metallicity variations; additional 
contributions may come from processes just beginning to become effective for ZAMS stars such as 
rotationally-induced mixing and planetary/circumstellar accretion. In older stars such as M67, 
the increase in scatter (and overall Li depletion) is then a product of processes efficiently acting 
on the main-sequence proper such as rotationally-induced mixing and/or photospheric accretion. 

Distinguishing specific mechanisms and their relative importance for Li depletion and scatter 
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at a given age will require continuing observational and theoretical efforts. Important advances on 
the theoretical front are at least three-fold: continued investigation of the role of magnetic fields in 
PMS Li depletion, realistic model atmospheres which include chromospheres, and detailed NLTE 
abundance calculations which employ these to extend extant sophisticated modeling attempts 
(e.g., Stuik et al. 1997). On the observational front, continued observations of Li in a variety of 
clusters spanning a range in age and metallicity are needed. We believe that particularly important 
observational work to be accomplished includes the determination of photometric periods in more 
cluster stars, detailed abundances of numerous elements (in particular, using both ionization 
sensitive and ionization insensitive features and elements) in cluster stars, quantification of even 
small "metal" (not just Fe) abundance spreads in cluster stars, and the association between 
planetary systems and parent star Li and light metal abundances. 

AK acknowledges support from NASA grant H-04630D to the University of Colorado. 
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TABLE 1 

Pleiades sikgle star sample 



Name 
(H II) 


V 


B-V 


V-I 


(S93a) 


(B-V) 


EW 
(B-V) 


Li 

(B-V) 


(V-I) 


EW 
(V-I) 


Li 

(V-I) 


<t(EW) 


<T(Li) 

(B-V) 


^(Li) 
(V-I) 


^(ALi) 
(B-V) 


fr(ALi) 
(V-I) 


Prot 
(d) 


25 


9.35 


0.44 


0.52 


6560 


6564 


57 


2.99 


6471 


57 


2.92 


5 


0.07 


0.07 


0.09 


0.09 




34 


11.84 


0.88 


0.92 


4900 


4928 


134 


1.94 


4996 


135 


2.03 


18 


0.14 


0.14 


0.11 


0.10 


6.553: 


152 


10.60 


0.645 


0.685 


5700 


5715 


150 


2.92 


57G0 


150 


2.97 


12 


0.10 


0.10 


0.08 


0.08 


4.12 ^ 


164 


9.40 


0.46 


0.54 


6480 


6474 


59 


2.95 


G37G 


58 


2.86 




0.07 


0.08 


0.09 


0.09 




193 


11.18 


0.7C 


0.80 


5320 


5339 


148 


2.51 


5352 


148 


2.52 


12 


0.11 


0.11 


0.07 


0.07 




250 


10.59 


0.645 


0.68 


5740 


5715 


141 


2.87 


5779 


141 


2.93 


12 


0.10 


0.10 


0.08 


0.08 




253 


10.59 


0.64 
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Fig. 1. — The color magnitude diagrams of our final non-binary Pleiades Li sample (open circles) 
and stars rejected as binaries (filled triangles). The Pleiades data is plotted assuming a distance 
modulus of 5.63 and with a 100 Myr isochrone from Pinsonneault et al. (1998). 
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Fig. 2. — The top panel shows the difference between our (B — V)- and (V — /)-based T^fi estimates 
versus Ha ffux ratio (relative to the stellar bolometric ffux) from Soderblom et al. (1993). The 
bottom panel shows the temperature difference versus the residual Ha flux ratio, which is the flux 
ratio less a fitted color-dependency. 
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Fig. 3. — Li abundances (with derived errors) vs. Tgfi from our {B — V) measures (top) and {V — I) 
measures. The weh-known dechning trend of Li abundance with decreasing T^s is fit with a fourth 
order Legendre polynomial (dashed line). 
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Fig. 4. — Same as Figure 3 except the differential (detrended) Li abundances and related errors 
are shown. Temperature bins used in considering the scatter of Li abundances as a function of T^s 
are labeled in both plots. 
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Fig. 5. Projected rotational velocity vsini versus {V — /)-bascd Li results for (a) all Tcff and 
absolute Li abundance (b) all Tgfj and differential detrended Li abundance (c) T^s = 4500 — 5500 
K and absolute Li abundance (d) T^s = 4500 — 5500 K and differential Li abundance. 
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Fig. 6. — Same as Figure 5 except rotational period, Prot, in days is plotted instead of vsmi. 
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Fig. 7. — Our (y — /)-bascd Li abundances versus stellar mass. The symbol size corresponds 
to rotational period. The largest circles represent stars with Prot < Id, the medium sized circles 
represent those with Pj-ot between Id and 5d, and the smaller circles indicate stars with Prot between 
5d and lOd. The dashed lines denote the mass range M = 0.7 — 0.9 Mq. 
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Fig. 8. — Our {V — /)-based differential Li abundances are plotted versus the Ca II infrared triplet 
flux ratio (top panel) and residual flux ratio (bottom), which is the flux ratio less a fltted color- 
dependency. 
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Fig. 9. — Differential A7699 K I equivalent widths ([observed— fitted] /fitted ) versus our (B — V)- 
based (top panel) and {V — /)-based differential Li abundances (bottom panel). 
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Fig. 10. — The {V — /)-based differential K I line strengths are plotted versus the Ca II infrared 
triplet flux ratio (top panel) and residual flux ratio (bottom panel) . 
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Fig. 11. — Differential A6717 Ca I equivalent widths versus our [B — l^)-based (top panel) and 
{V — /)-based (bottom panel) differential Li abundances. 



